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ABSTRACT 

The Kepler mission has yielded the discovery of eight circumbinary systems, all found around 
eclipsing binaries with periods greater than 7 d. This is longer than the typical eclipsing binary 
period found by Kepler, and hence there is a dearth of planets around the closest binaries. In 
this paper we suggest how this dearth may be explained by the presence of a distant stellar ter¬ 
tiary companion, which shrunk the inner binary orbit by the process of Kozai cycles and tidal 
friction, a mechanism that has been implicated for producing most binaries with periods be¬ 
low 7 d. We show that the geometry and orbital dynamics of these evolving triple-star systems 
are highly restrictive for a circumbinary planet, which is subject itself to Kozai modulation, on 
one hand, and can shield the two inner stars from their Kozai cycle and subsequent shrinking, 
on the other hand. Only small planets on wide and inclined orbits may form, survive and allow 
for the inner binary shrinkage. Those are difficult to detect. 

Key words: binaries; close, eclipsing - astrometry and celestial mechanics: celestial mechan¬ 
ics - planets and satellites: detection, dynamical evolution and stability - methods: analytical, 
numerical 


1 INTRODUCTION 

The first two decades of exoplanetary science have yielded many 
surprising results. Not only do most stars host orbiting planets (Cas- 
san et al. 2012), but planets are often found in unexpected locations 
and with unexpected properties. For example, hot-Jupiters continue 
to pose significant theoretical challenges (Triaud et al. 2010; Mad- 
husudhan et al. 2014), while Super-Earths were predicted not to 
form, and yet they are some of the most abundant planets known to¬ 
day (Howard et al. 2010; Mayor et al. 2011). Planets also have been 
found in binary star systems orbiting one (e.g., 16 Cygni, Cochran 
et al. 1997) and two (e.g., Kepler-16, Doyle et al. 2011) stars. The 
parameter space of non-discoveries is shrinking fast. The field has 
therefore evolved to a state where an absence of planets is just as 
telling as a new discovery. 

One conspicuous absence is seen in the Kepler circumbinary 
planets (CBPs). So far there have been ten transiting CBPs discov¬ 
ered by Kepler orbiting eight eclipsing binaries (EBs), including the 
three-planet system Kepler-47 (Orosz et al. 2012b). It was pointed 
out by Welsh et al. (2014a) that all of the planets have been found 
orbiting EBs of periods between 7.4 and 40 d, despite the median 
of the EB catalog being 2.7 d. The discoveries have therefore been 
made on the tail of the EB period distribution. The analyses of Arm- 
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strong et al. (2014) and Martin & Triaud (2014) suggest that this 
dearth of planets is statistically significant. 

In this paper we show how this absence of planets may be a 
natural consequence of close binary formation, formalising an ex¬ 
planation that was briefly proposed by Welsh et al. (2014a), Arm¬ 
strong et al. (2014), Martin & Triaud (2014) and Winn & Fabrycky 
(2014). A popular theory proposes that most very close binaries 
(< 7 d) are initially formed at wider separations. The binaries sub¬ 
sequently shrink under the influence of a misaligned tertiary star 
and a process known as Kozai cycles with tidal friction (KCTF) 
(e.g., Mazeh & Shaham 1979; Fabrycky & Tremaine 2007). The 
dynamics of such evolving triple star systems impose strict stabil¬ 
ity constraints on any putative circumbinary planetary orbits around 
the inner binary. We will demonstrate that these constraints act to 
either completely inhibit planetary formation and survival or re¬ 
strict it to small planets on wide, misaligned orbits. 

We point out that Munoz & Lai (2015) and Hamers et al. 
(2015b) have conducted independent analyses with complementary 
techniques, reaching the same general conclusions as this paper. 

The plan of this paper is as follows: In Sect. 2 we discuss 
some of the relevant trends seen in the Kepler circumbinary plan¬ 
ets. Sect. 3 is a review of the theory and observations related to 
close binary formation via KCTF. In Sect. 4 we summarise some 
important theoretical aspects of three-body dynamics and stability 
in binary star systems, and then extend the discussion to the triple 
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Table 1. The Kepler transiting circumbinary planets 


Name 

Ml 

(Mol 

M2 

(Mol 

^in 

[AU] 

Pm 

Id] 

^in 

Pp 

[P®] 

<3p 

[AU] 

[d] 

ep 

A/p,m 

[deg] 

^crit 

[AU] 

Reference 

16 

0.69 

0.20 

0.22 

40.1 

0.16 

8.27 

0.71 

228.8 

0.01 

0.31 

0.64 

Doyle etal. (2011) 

34 

1.05 

1.02 

0.23 

28.0 

0.52 

8.38 

1.09 

288.8 

0.18 

1.86 

0.84 

Welsh et al. (2012) 

35 

0.89 

0.81 

0.18 

20.7 

0.14 

7.99 

0.60 

131.4 

0.04 

1.07 

0.50 

Welsh et al. (2012) 

38 

0.95 

0.26 

0.15 

18.8 

0.10 

4.35 

0.47 

106.0 

0.07 

0.18 

0.39 

Orosz et al. (2012a) 

47b 

1.04 

0.36 

0.08 

7.4 

0.02 

2.98 

0.30 

49.5 

0.04 

0.27 

0.20 

Orosz et al. (2012b) 

47d 

1.04 

0.36 

0.08 

7.4 

0.02 

— 

0.72 

187.3 

— 

— 

0.20 

Orosz (in prep) 

47c 

1.04 

0.36 

0.08 

7.4 

0.02 

4.61 

0.99 

303.1 

<0.41 

1.16 

0.20 

Orosz et al. (2012b) 

PH-1/64 

1.50 

0.40 

0.18 

20.0 

0.21 

6.18 

0.65 

138.5 

0.07 

2.81 

0.54 

Schwamb et al. (2013); Kostov et al. (2013) 

413 

0.82 

0.54 

0.10 

10.1 

0.04 

4.34 

0.36 

66.3 

0.12 

4.02 

0.26 

Kostov et al. (2014) 

3151 

0.93 

0.19 

0.18 

27.3 

0.05 

6.17 

0.79 

240.5 

0.04 

2.90 

0.44 

Welsh et al. (2014b) 


star case. In Sect. 5 we run n-body simulations on a set of example 
systems to test planet stability and binary shrinkage. We then dis¬ 
cuss the protoplanetary disc environment and the effects of planet 
migration and multi-planetary systems in Sect. 6. In Sect. 7 we con¬ 
clude by summarising the general argument and looking ahead to 
future complementary observations. 


2 CIRCUMBINARY PLANETS DISCOVERED BY Kepler 

The Kepler mission provided four years of near-continuous pho¬ 
tometry of roughly 200,000 stars. As of May 8, 2015 there have 
been 2773 Kepler eclipsing binaries identified (Prsa et al. 2011; 
Slawson et al. 2011, Kirk et al. in prep)'. A search for additional 
transit signals has led to the discovery of eight circumbinary sys¬ 
tems, one of which contains three planets (Kepler-47) (Welsh et al. 
2014a and references there-in). 

In Table 1 we summarise the parameters of the 10 transiting 
CBPs discovered by Kepler. We list the masses of the two stellar 
components. Mi and M 2 , the semi-major axis, ai,,, period, Pi„, and 
eccentricity, of the inner binary orbit. We then list the planetary 
parameters: its radius, R^, semi-major axis, flp, period, Pp, eccen¬ 
tricity, gp and mutual inclination with respect to the inner binary 
orbital plane, A/pi„. We then list the inner-most stable orbit for 
a circumbinary planet, flcricCB, calculated according to Holman & 
Wiegert (1999) (see Sect. 4.1.1 for details). The last column is the 
discovery paper reference, including the two independent discover¬ 
ies of PH-l/Kepler-64 by Schwamb et al. (2013)^ and Kostov et al. 
(2013). The planet masses are not listed because they are generally 
poorly constrained. 

In the small sample of CBPs some preliminary trends have 
been identified: 

(i) The planets have all been discovered orbiting eclipsing bi¬ 
naries of relatively long periods. In Fig. la we demonstrate this 
feature by plotting a histogram of the Kepler eclipsing binary peri¬ 
ods, which has a median of 2.7 d (blue solid vertical line). The red 
dashed vertical lines indicate the binary periods around which plan¬ 
ets have been found, of which the shortest is = 1 .A d (Kepler- 
47). Armstrong et al. (2014) analysed the Kepler lightcurves using 
an automated search algorithm and debiasing process, and came to 

' The latest version of the catalog can be found online at 
http://keplerebs.villanova.edu, maintained by Prsa et al. at Villanova 
University. 

^ The planet hunters consortium at http://www.planethunters.org/. 


the conclusion that planets are significantly rarer around binaries of 
periods between 5 and 10 d than around wider binaries. Martin & 
Triaud (2014) created synthetic circumbinary distributions and sim¬ 
ulated eclipses and transits observable by Kepler. They concluded 
that the number of discoveries would have been roughly doubled 
if the circumbinary abundance did not drop around binaries with 
periods shorter than 5 d. 

(ii) Most of the circumbinary planets are found near the dynam¬ 
ical stability limit, which we demonstrate in Fig. lb. Welsh et al. 
(2014a) proposed that this was either the result of a physical pile- 
up of planets or detection biases, although Martin & Triaud (2014) 
argued against the latter explanation. 

(iii) The systems are all close to coplanarity. Out of the ten cir¬ 
cumbinary planets found so far, the most misaligned is Kepler-413 
with A/pin = 4° (Kostov et al. 2014). The mean misalignment is 
1.7°. This is however an observational bias imposed by the require¬ 
ment of consecutive transits in the detection method. Significantly 
misaligned systems (> 10°) would either avoid transiting or only 
do so irregularly (Martin & Triaud 2014), hindering the ability to 
detect them. 

(iv) The planets have only been found within a narrow size range 
between 3 and 8.3 R^,. The lack of Jupiter-size and larger planets 
might be due to the relative rarity of such big planets and the small 
number of detections. The paucity of small planets is probably a 
result of the large variations in circumbinary transit timing and du¬ 
ration, making it difficult to phase-fold the photometric data and 
identify shallow transits. 

(v) Preliminary calculations indicate that the circumbinary 
abundance, within the present realm of detections, is ~ 10% (Arm¬ 
strong et al. 2014; Martin & Triaud 2014), which is similar to the 
13% abundance of Neptune-mass and heavier planets around sin¬ 
gle stars derived from radial velocity surveys (Mayor et al. 2011). 
This circumbinary abundance was calculated based on only the 
near-coplanar population that Kepler is sensitive to. If there exists 
a presently-hidden population of misaligned circumbinary planets 
then this would increase their overall abundance. 


3 CLOSE BINARY FORMATION WITH A TERTIARY 
STELLAR COMPANION 

A common belief is that short period binaries are not formed in 
situ, because the fragmenting protostellar cloud contains an excess 
of angular momentum with respect to the orbital angular momen¬ 
tum of a very close binary (e.g.. Bate 2012). These binaries were 
likely formed much farther apart and subsequently evolved to their 
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Figure 1. Circumbinaiy planets discovered so far via Kepler transits: (a) A histogram of periods of the Kepler EB catalog, with a blue solid vertical line 
indicating the median of this population. Red dashed vertical lines indicating the binary periods known to host circumbinary planets, (b) Planet semi-major 
axis vs critical semi-major axis (Eq. 4), where the blue shaded region corresponds to unstable orbits. Systems have been labeled by their Kepler number. The 
outer planets of Kepler-47 are presented by square symbols. 


(a) 



(b) 



Figure 2. Example of Kozai cycles of (a) and A/in,out (b) for Pin = 100 
d, Pout = 338 yr, Mi = M 2 = M 3 = Mq and A/in out,init = 67°. Both the 
inner and outer binaries are on initially circular orbits. 


present state. An alternate theory for the formation of the closest of 
binaries (< 7 d) is the combination of tidal interactions and an in¬ 
clined perturbing tertiary star. This process is canonically known as 
Kozai cycles with tidal friction (KCTF) and was first proposed by 
Mazeh & Shaham (1979), and later studied in detail by Eggleton & 
Kisseleva-Eggleton (2006), Eabrycky & Tremaine (2007) and Naoz 
& Eabrycky (2014) among others. In this section we summarise this 
shrinking mechanism. 

3.1 The Kozai modulation 

To a first approximation, a hierarchical triple star system can be 
modelled as an inner binary of two stars and an outer binary com¬ 
posed of the inner binary, located at its centre of mass, and the outer 
tertiary star. Both binaries move on Keplerian orbits which we de¬ 
fine using osculating orbital elements for the period, P, semi-major 
axis, a, eccentricity, e, argument of periapse, a>, and longitude of 
the ascending node, Q., where we denote the inner and outer bina¬ 
ries with subscripts “in” and “out”, respectively. The two orbits are 
inclined with respect to each other by A/in out- The inner binary stars 
have masses Mj and M 2 and the tertiary mass is M 3 . 

In the quadrupole approximation of the Hamiltonian of the 
system the tertiary star remains on a static orbit whilst its pertur¬ 
bations induce a nodal and apsidal precession on the inner binary 
(variations in fljn and Wjn, respectively). If the two orbital planes 
are initially misaligned by more than a critical value'^. A/in,out,init > 
39.2°, there is a variation of gin and A/in,out on the same period as 
the precession of Win, and a nodal precession of flin on twice this 

^ This critical mutual inclination is for an initially circular inner binary. 
The critical value decreases with increasing initial eccentricity in the inner 
binary. 
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Figure 3. Geometi'y of a close binary orbited by a planet and a distant mis¬ 
aligned stellar companion. The orbits are drawn with respect to the centre of 
mass of the inner binary. All orbits are actually circular but misaligned with 
respect to our line of sight. The planet is coplanar with the inner binary. 


period (see Takeda et al. 2008 for a more in-depth summary). The 
quadrupole approximation works best when aout/flin is very large, 
and we will be considering systems like this. 

An initially circular inner binary reaches an eccentricity of 


amplitude of the tertiary’s perturbations decreases and the Kozai 
effect is suppressed. 

If the inner binary stars are close enough then apsidal preces¬ 
sion due to general relativity and tidal and rotational bulges can 
suppress the Kozai modulation (Wu & Murray 2003; Fabrycky & 
Tremaine 2007). Alternatively, a sufficiently close and massive cir- 
cumbinary planet may also suppress the Kozai modulation by in¬ 
ducing a competing apsidal advance, an effect we consider later in 
Sect. 4.2.3 and has been studied by Hamers et al. (2015a). 


3.3 Tidal shrinkage 

Assuming the Kozai cycles have not been suppressed, the high ec¬ 
centricity excursions lead to close periapse passages between the 
two stars of the close binary. During these close encounters there 
may be a strong tidal interaction, which can dissipate orbital energy 
that causes the inner binary orbit to both shrink and circularise. The 
final semi-major axis can be approximated as twice the periapse 
distance at closest approach. 


^in,final “ 2flin,init(l ^in.max). 


(3) 


^in.max 


1 2 ij,n 


1/2 


( 1 ) 


whilst the mutual inclination varies between its initial value and the 
critical value 39.2°. These variations together constitute the Kozai 
cycles (Lidov 1962; Kozai 1962)^'. The timescale of Kozai cycles is 


^Kozai,in ~ 


2 Pi,,, Ml+M2 + 


37t P„ 




(1 


1 


3/2 


( 2 ) 


(e.g., Mazeh & Shaham 1979; Fabrycky & Tremaine 2007). The 
period of Kozai cycles equals TKozai.in multiplied by a factor of order 
unity (Ford et al. 2000). 

A numerical example is given in Fig. 2, demonstrating that 
can be excited to very high values. This simulation, and all other 
simulations in this paper, was created using the n-body code RE¬ 
BOUND (Rein & Liu 2012)^. The code uses a 14/15th-order inte¬ 
grator which is not inherently symplectic but uses an adaptive time 
step to minimise error propagation generally down to machine pre¬ 
cision. Consequently, it preserves the symplectic nature of Hamil¬ 
tonian systems better than typical symplectic integrators (Rein & 
Spiegel 2015). The code is also suitable for high-eccentricity or¬ 
bits, like those experienced during Kozai cycles. All parameters for 
the integration were kept at default. 


3.2 Suppression of Kozai cycles 

Kozai cycles are the result of small perturbations induced by the 
tertiary star which build up coherently over the apsidal precession 
period. If there is an additional secular perturbation causing an ap¬ 
sidal precession of the binary on a shorter timescale, then the co¬ 
herent eccentricity modulation is partially lost. Consequently, the 


although this formula is only valid for when the two stars are close 
enough for tidal forces to be significant (Ford & Rasio 2006). 

Observational studies of binary and triple star systems (e.g., 
Tokovinin 1993, 2004; Tokovinin et al. 2006; Tokovinin 2008, 
2014a,b) are consistent with this shrinkage scenario. Tokovinin et 
al. (2006) discovered that 96% of binaries with less than 3 d peri¬ 
ods are surrounded by a tertiary stellar companion. This abundance 
drops to 34% for binaries with periods longer than 12 d. 


4 PLANETARY DYNAMICS IN MULTI-STELLAR 
SYSTEMS 

There are two possible planetary orbits in binary star systems: 

(i) a circumbinary orbit around both stars or 

(ii) a circumprimary orbit around one of the two stars^. 

These orbits are also sometimes referred to as p-type and s-type 

orbits, respectively, but we will not use this nomenclature. 

There is an even greater variety of planetary orbits in triple star 
systems (see Venier & Evans 2007 for a classification scheme). 
In this paper we are only concerned with the scenario illustrated 
in Eig. 3, where the planet orbits the inner binary and the third 
star is on the periphery. The planet may be considered as both on 
a circumbinary orbit around the inner binary, and on a circumpri¬ 
mary orbit with respect to the outer binary, where the planet orbits 
a “star” that is composed of the inner binary. In Sect. 4.1 we anal¬ 
yse both types of orbits in the context of binary star systems, before 
combining them in the context of triple stellar systems in Sect. 4.2. 


This effect was first published (in Russian) by Lidov but was indepen¬ 
dently discovered by Kozai later that year. It should arguably be called a 
Lidov-Kozai cycle but it is more commonly refeired to as simply a Kozai 
cycle, a convention which we follow in this paper. 

^ This easy to use code can be downloaded freely at 
http://github.com/hannorein/rebound. 


® Technically a circumprimary orbit only refers to when a planet orbits 
the bigger of the two stars and an orbit around the smaller star is a cir- 
cumsecondary orbit, but for simplicity in this paper we will use the term 
circumprimary to refer to either case. 
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4.1 Planets in binary star systems 

4.1.1 Stability 

There are restrictions on where a planet may orbit stably in a binary 
star system, primarily as a function of the binary’s semi-major axis 
(e.g., Dvorak 1986; Holman & Wiegert 1999). Planets orbiting on 
the wrong side of the stability limit are generally ejected from the 
system by a process of resonance overlap (Mudryk & Wu 2006). 
Assuming circular circumbinary orbits Holman & Wiegert (1999) 
used an empirical fit to n-body simulations to define the stability 
limit using a critical semi-major axis. 


rtcrit.CB “ r/in(1.60 + fi.lOCin +4.12^in 4.27CiiiPin 

- 2.224-5.09;rf„ + 4.61c?XX (4) 

where = MiKMi + M 2 ) is the mass ratio. We have excluded the 
fit uncertainties that the authors included in their equation, which 
are of the order 5%. Planets must orbit beyond this critical semi¬ 
major axis, lest they are ejected from the system. 

In the alternative case of a circular circumprimaty orbit Hol¬ 
man & Wiegert (1999) similarly calculated a stability criterion of 



-30 -20 -10 0 10 20 30 

X [AU] 


^crit.CP rtoij[(0.46 0.63Cout "H O.fiOCoutPout 

-I- 0.154t - 0.204tPout), (5) 

whereyUout = A/ 3 /(Mi-(-M 2-1 M 3 ) is the mass ratio. Equation 5 does 
not account for eccentric planets. In this configuration the planet 
must orbit within this critical semi-major axis in order to maintain 
stability. Both Eqs. 4 and 5 were derived for a binary eccentricity 
between 0 and 0.7-0.8, and a binary mass ratio between 0.1 and 
0.9. 

The analysis of Holman & Wiegert (1999) was restricted to 
coplanar orbits. Doolin & Blundell (2011) showed that in the cir¬ 
cumbinary case flcrit is only a weakly decreasing function with mu¬ 
tual inclination, and that stable circumbinary orbits are possible at 
all misalignments, including retrograde orbits. 


4.1.2 Secular evolution 

A circumbinary planet experiences a nodal and apsidal precession, 
of tUp and Dp, respectively, both of which occur at approximately 
the same rate but in opposite directions (e.g., Lee & Peale 2006). 
The timescale of this precession is 


('■'1 

(Ml -t-Mz)^ 1 


l^tuj 

1 Ml M 2 COsA/p3n 


( 6 ) 


which is adapted from Earago & Laskar (2010) to be in terms 
of orbital periods. Eor circular binaries the mutual inclination A/ 
remains constant, but for eccentric binaries it may vary slightly 
(Earago & Laskar 2010; Doolin & Blundell 2011). 

A circumprimary orbit will also experience a precession tUp 
and Dp. If the orbit is misaligned by at least 39.2° with respect to 
the outer binary it will undergo Kozai cycles, as was seen for stellar 
binaries in Sect. 3.1, with a timescale 


TKo. 


2 P4 Ml -t- Mp -h M 2 


2>n P„ 


M 2 


(1 - ^L) 


3/2 


(7) 


Note that Kozai cycles only apply to circumprimary orbits, not to 
circumbinary ones. 


Figure 4. Stability limits in a triple star system from Eqs. 4 and 5, for 
Ml = M 2 = M 2 , = Mq, Qin = 0.5 AU, flout = 30 AU, = Sout = 0 and all 
orbits are coplanar. The two blue shaded regions are unstable for planetary 
orbits and the white region in between is stable. Note that this does not take 
into account eccentric planets. 


4.2 Planets in triple star systems 

4.2.1 Stability 

Because a planetary orbit in a triple star system has both circumbi¬ 
nary and circumprimary characteristics, it must obey both stability 
constraints Qcrit.cB (Eq. 4) and flcrit.cp (Eq. 5). If flout/oin is suffi¬ 
ciently large then flcrit.cp > 'tcrit.cB. and hence there is a region of 
stable planetary orbits. In Fig. 4 we illustrate these combined sta¬ 
bility limits, using a simple example of ain = 0.5 AU, flout = 30 AU, 
equal mass stars and coplanar and circular orbits. 

Verrier & Evans (2007) analysed this scenario with copla¬ 
nar and circular orbits and concluded that this joint application of 
Eqs. 4 and 5 is a valid description of planetary stability in triple 
star systems. This implies that the perturbations from the inner and 
outer binaries act somewhat independently. This assumption was 
seen to become less valid when the inner and outer binary orbits 
were made eccentric, resulting in a smaller stability region than 
what is defined by Eq. 4 and Eq. 5. This additional instability was 
attributed to combined perturbations from the ensemble of three 
stars. Verrier & Evans (2007) did not test planet eccentricity, which 
would no doubt make the stability constraints even more restrictive. 


4.2.2 Secular evolution 

A planet orbiting within the stability region of a misaligned triple 
system is perturbed by two competing secular effects: precession 
due to the inner binary and Kozai cycles due to the outer binary. 
Kozai cycles in the planet may be suppressed by the precession in¬ 
duced by the inner binary, like what was discussed in the case of 
triple star systems in Sect. 3.2. An approximate criterion for sup¬ 
pression of Kozai cycles is 
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Figure 5. Competing secular timescales on the planet: circumbinary pre¬ 
cession (Eq. 6), induced by the inner binary on the planet, and Kozai cycles 
(Eq. 7), induced by the third star at 100 AU. The circumbinary precession 
timescale is shown in different shades of blue for different values of Pi„, 
starting from 100 d (dark blue) to 5 d (light blue). The Kozai timescale is 
shown for Cout = 0 (maroon) and 0.25 (red) and is independent of P,„. This 
plot covers the simulated range of Up between 1.5 and 10 AU. 


^prec.p ^Kozai.p- (8) 

In reality there is likely to be a transition between the two com¬ 
peting secular effects when their timescales are similar. Verrier & 
Evans (2009) showed that the inner binary can “protect” the planet 
through this suppression of Kozai cycles. On the other hand, planets 
left to undergo high-amplitude Kozai cycles were seen by Verrier & 
Evans (2009) to become unstable rapidly. Some alternative means 
of suppressing Kozai cycles, such as general relativistic precession 
and tidal and rotational bulges in the stars, are not applicable here. 
This is because these effects are only non-negligible for planets 
very close to the stars, and such orbits are prohibitively unstable 
for circumbinary planets. The Kozai effect may also be suppressed 
by perturbations from an additional planet, but we do not consider 
multi-planetary systems in this analysis (see Sect. 6.3 for a short 
discussion on the effects of additional planets). 

We have not followed the analytic derivation to the octopole 
order, which can lead to secular resonance between different pre¬ 
cession terms resulting in ~ 0.1 eccentricity fluctuations over semi¬ 
major axis ranges of ~ 10% (Naoz et al. 2013). This effect may lead 
to additional unstable regions in the CBP parameter space. 

As the binary orbit shrinks under KCTE the planet precession 
timescale increases as Tprec.p oc whilst Txozai.p remains con¬ 

stant. This means that the relative influence of Kozai cycles on the 
planet will increase and may start to dominate. 


4.2.3 KCTF suppression due to planetary mass 

Just as Kozai cycles of the planet may be suppressed by the apsidal 
motion induced by the host binary, the inner binary might fail to 
excite to high eccentricity as a result of apsidal precession due to 
the mass of the planet. In effect, if the planet’s tidal force exceeds 
that of the tertiary star, the host binary processes too quickly for 
the tertiary to secularly excite its eccentricity. This is possible even 
for small planetary masses, because the planet is much closer to 
the binary than the third star is, and tidal effects scale as separation 
cubed. 


For Kozai cycles analysed to quadrupole order, we follow Fab¬ 
rycky & Tremaine (2007) Sect. 3.2 to compute the maximum ec¬ 
centricity of the binary in the presence of additional precession. For 
a triple system, there are two constants of motion: 

F' = -2 - 3e^ -H (3 -t 12^2 - 15e cos^ w) sin^ AI , (9) 

H' = (1-e^)‘/2cosA/, (10) 


where the former is a non-dimensionalised form of the Hamilto¬ 
nian term defining the interaction between the two orbits, and the 
latter is a non-dimensionalised form of the z-component of angular 
momentum^, which is the component of the angular momentum of 
the binary along the total angular momentum of the system. In this 
context, e and oj are for inner binary and AI is the mutual inclina¬ 
tion of the two orbits. 

In the case of a quadruple system, we may consider each level 
of the hierarchy as an orbit of a body around the centre of mass of 
the interior bodies. Averaging over the orbital trajectories gives the 
following averaged Hamiltonian (e.g.. Ford et al. 2000): 


{T) 

= 

*^11 "I" "1" "I" "1" {*^, 001 ) "1" (*^ 11 , 001 ); 

(11) 

< 9 t„> 


GA/| A/2 

(12) 

(Tp) 


G(Mi -H M2)Mp 

( 13 ) 

<r„ut> 


G(Mi +M 2 + Mp)M3 

2^Zout 

( 14 ) 

<? 7 n.p> 


GM1M2MP al 

M 1 +M 2 SalJl-el,fF 

( 15 ) 


X 

(2 - 1 - 3el -(3 + I2el - 15 ef„ cos^ win.p) sin^ A/,„,p), 

(Tp.out) 

= 

G(Mi -1- M 2 )MpM 3 al 

(M, + M2+ Mp 8ap(l - elfF 

( 16 ) 


X 

(2 - 1 - 3 ep - (3 -b 12ep - ISCp cos^ tUp) sin^ A/p,o„t^ 


(*^n,out) 

= 

GM1M2M3 al 

M,+M2 Sal,d-el,fF 

( 17 ) 


X 

(2 - 1 - 3el - (3 -b 12 ef„ - 15 ef„ cos^ ojin) sin^ A/i„.„ut). 


Suppose the binary begins on a circular orbit, with a copla- 
nar exterior planet (A/pi„ = 0). The planet’s orbit continues on a 
nearly coplanar orbit, for the same reason that its Kozai cycles are 
suppressed, as just described in Sect. 4.2.2. We no longer have con¬ 
servation of F', but rather conservation of ((F). In the secular prob¬ 
lem, without tidal dissipation, a^, Op, and flom are all conserved, 
so the Keplerian Hamiltonian terms, Eqs. 12 -14, are individually 
conserved. The other sub-components can trade energy. 

Now define 


"p «L(i 






M3 aid-elf F 
(Ml -H Mzj^Mp 

MiM2(Mi -{- M 2 -t- Mp) \a 


(18) 

(19) 


which measure the numerical dominance of the different terms. For 
typical circumbinary planets,/? « 1 because Mp < M^ia^laiff. In 
that case, ((Tj,,out) is negligible compared with ((Fn.p) and ((Fn,out)> 
and we can find a relatively simple closed-form solution of the ec¬ 
centricity maximum of the binary. This is accomplished by making 


^ This corrects a misprint in equation 17 of Fabrycky & Tremaine (2007). 
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flp = 1.5 AU 


Qp = 2 AU 



flp = 9 AU 


Inner binary 



Figure 6. First seven figures: evolution over 100 Myr of the planet eccentricity in a primordial triple star system, where din = 0.4827 AU (Pin = 100 d), 
flout = 100 AU (Pout = 2.52 X 10^ d) and ^out = 0 (dark blue) and 0.25 (light blue) and the inner binary starts on a circular orbit. The planet is massless and 
starting on an orbit that is circular and coplanar with the inner binary. Each subplot shows a different value of dp (see Table. 2 for all the simulation parameters). 
Bottom right figure: evolution of the inner binary eccentricity for ^out = 0 (dark blue) and 0.25 (light blue). 
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Figure 7. First seven figures: evolution over 100 Myr of the mutual inclination between the planet and the inner binai-y in a primordial triple star system, where 
Qin = 0.4827 AU (Pin = 100 d), Uout = 100 AU (Pout = 2.52 x 10^ d) and eout = 0 (dark blue) and 0.25 (light blue) and the inner binary starts on a circular 
orbit. The planet is massless and starting on an orbit that is circular and coplanar with the inner binary. Each subplot shows a different value of Up (see Table. 2 
for all the simulation parameters). Note the changing scale on the y-axis as a function for different Up. Bottom right figure: evolution of the mutual inclination 
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the substitution 


sin A/in,out “ 1 cos A/in,out “ 1 cos A/in,out.init(l ^in) ; (20) 


due the initial condition ein.init = 0 and the conservation of H'. The 
result is an implicit equation of Cin versus Wj,,. The maximum ec¬ 
centricity occurs at Wi„ = ±nl2, which we insert and manipulate to 
yield: 


^in,max 


3-0- 


- cos'^ A/i„ 


1/2 


( 21 ) 


This expression recovers Eq. 1 in the limit o —> 0, and shows that 
the maximum value of eccentricity cycles declines with increasing 
planet mass until no eccentricity is excited for masses above: 


Mp,cri, = 3M3 




( 22 ) 


Alternately, for known CBP systems, we may view this as a 
limit on the properties of a tertiary that can cause the binary to 
shrink by Kozai cycles with tidal friction: 


I _ Ml _\ =_ Me _ (23) 

where the left and right hand sides quantify the strength of the tidal 
force from the tertiary star and planet, respectively. Kozai is sup¬ 
pressed on the inner binary when the left hand side is larger. The 
only known CBP that has a third distant companion is PH-l/Kepler- 
64; this companion is itself a binary (Ba-Bb) of total mass ~ 1 .5Mq. 
The planet’s mass is not known, but its radius of ~ 6Rq suggests its 
mass is at least of order Neptune’s, i.e. 154/®. The critical value 
of the planet’s tidal parameter, the right hand side of the last equa¬ 
tion, is therefore 6 x lO^^MoAU^^. On the other hand, the projected 
separation of Ba-Bb to the host binary is 1000 AU. Suppos¬ 
ing its eccentricity is very high and so its semi-major axis is only 
flout = 500 AU, then the critical value of the parameters is reached 
only for gout = 0.998. This is too high to allow the planet to remain 
dynamically stable. Hence we essentially rule out the possibility 
that this companion will excite Kozai cycles in the planet-hosting 
binary, because the planet processes its host too fast to allow this to 
happen. 


5 N-BODY SIMULATIONS OF A SET OF EXAMPLE 
SYSTEMS 

In this section we use n-body simulations to determine regions 
where planets may be able to form and survive in evolving triple 
star systems, and whether or not the binary will be allowed to shrink 
under KCTF. We take a stable orbit to be a necessary but probably 
not sufficient condition for planetary formation. We are therefore 
considering the most optimistic scenario possible. Additional ef¬ 
fects are considered afterwards in Sect. 6 . The entire argument is 
summarised in Sect. 7.1. 


5.1 A primordial inner binary 

For our first suite of n-body simulations we constructed an example 
primordial triple star system, with an interior circumbinary planet, 
using the parameters listed in Table 2. For the inner binary we 
tested a single set of parameters, where the orbit was arbitrarily 
chosen at 100 d and the stellar masses were the average values for 
the circumbinary systems discovered by Kepler {Mi = IM® and 


Table 2. Parameters for the primordial triple stai' system and circumbinary 
planet 


inner binary (one set of parameters) 


Ml 

1 M® 

M 2 

0.5 M® 

^in 

0.4827 AU 

Pm 

100 d 

^in 

0 

circumbinary planet (18 sets of parameters) 


Mp 

0 

£3in 

1.5 — lOAU 

Pm 

548 — 9430 d 

gp 

0 


0 

outer binary (two sets of parameters) 


M 3 

0.6 M® 

^out 

100 AU 

^out 

2.52 X 10^ d 

^out 

0, 0.25 

‘^An,out 

74° 


M 2 = O.SMq). The inner eccentricity was set to zero initially, al¬ 
though it rose dramatically during Kozai cycles. 

The tertiary star was given a mass M 3 = O.fiM®, correspond¬ 
ing to the median observed outer mass ratio ^out = 0-4, and placed 
at a typical separation of 100 AU (Tokovinin 2008). The mutual 
inclination A/n.out = 74° was chosen such that the inner binary 
would ultimately shrink to a 5 d period, according to Eq. 3, where 
^in.max = 0-93 (Eq. 1). We tested two different tertiary star orbits 
where we only changed gout = 0 and 0.25, to test the effect it has on 
planetary stability. We also tested e^^i = 0.5 but found that almost 
all systems were unstable in this configuration so we have omitted 
these results for simplicity. The requirement that gout < 0.5 for sta¬ 
bility is restrictive because observations show this to be a typical 
outer eccentricity (Tokovinin 2008). 

Orbiting around the inner binary we placed a massless planet 
on a circular and coplanar orbit, like in Fig. 3. The behaviour and 
stability of the planet are largely functions of the flp, as this de¬ 
termines the relative perturbing strengths from the inner and outer 
binaries. For the simulations we tested 18 different values of flp be¬ 
tween 1.5 and 10 AU, in steps of 0.5 AU. 

The range of semi-major axes was chosen based on the relative 
Kozai and precession timescales on the planet. We plot these com¬ 
peting timescales in Fig. 5 as a function of flp. The Kozai timescale 
in the planet is a monotonic decreasing function of flp (Eq. 7)*. 
We plot TKozai.p the two values of gout. On this logarithmic scale 
gout does not have a discernible affect on TKozai.p- The circumbinary 
precession timescale (Eq. 6 ), on the other hand, increases with flp. 
For Pin = 100 d there is a turnover between the two timescales at 
approximately = 1 AU. For farther out planetary orbits we ex¬ 
pect the planet to undergo Kozai cycles and obtain an eccentricity 
ep.max = 0.93, leading to ejection. For smaller flp, we expect Kozai 
to be suppressed and the planet to remain stable, as long as it is 
beyond the circumbinary stability limit in Eq. 4. 

Using the REBOUND code we ran n-body integrations for the 
total of 36 simulations. These were run over 100 Myr, which lets 


“ To apply this equation to a triple star system one replaces Mi by Mj -H M 2 
and M 2 by M 3 . 
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us cover almost 100 Kozai timescales (TKozai.in = 1.17x 10® yr from 
Eq. 2). When the tertiary star is eccentric there is also a longer 
octupole timescale, 


ttout(l ^OUt) 
^out^in 


1/2 


^ Kozai,i 


(24) 


(Antognini 2015). In our example Toct.in = 3.27 x lO’ yr, and hence 
we are covering a couple of these timescales, but may be missing 
some dynamical effects that develop over many timescales. In our 
simulations a planet was considered unstable if it were ejected from 
the system. 

In Fig. 6 we plot over time for a selection of simulations: 
flp = 1.5, 2, 3.5, 4, 7, 8 and 9 AU. In each plot we show the result 
for fiout = 0 (dark blue) and 0.25 (light blue). For reference, in the 
bottom right figure we show the inner binary over time. The max¬ 
imum fin is negligibly different in the two cases: fm.max = 0.9239 
for Com = 0 and ein max = 0.9241 for = 0.25 but the Kozai mod¬ 
ulation period changes significantly: ~ 10 yr for a circular tertiary 
and ~ 5.5 yr for an eccentric tertiary. 

For flp = 1.5 AU neither configuration is stable, as in each 
case the planet is ejected as soon as the inner binary reaches its 
eccentricity maximum for the first time. When the planet is a little 
farther from the inner binary, at 2 AU, and eout=o, its eccentricity 
rises to 0.1 periodically, coinciding with when the inner binary un¬ 
dergoes its Kozai modulation, but the planet nevertheless remains 
stable over 100 Myr. This roughly calculated inner stability limit 
coincides well with Eq. 4 from Holman & Wiegert (1999), where 
flcrit.cB = 1.94 AU for ein = ein.max = 0.94. When gout is increased 
to 0.25 the planet is no longer stable at 2 AU. In fact the shortest- 
period planet that was stable over 100 Myr for = 0.25 was 
flp = 3.5 AU, although even in this configuration the planet ob¬ 
tained significant eccentricity that may have led to instability on a 
longer timescale. For Op = 4 AU there is no large eccentricity varia¬ 
tion and the results gout = 0 and 0.25 are almost identical. This may 
be considered the inner stability limit for fiom = 0.25. 

The inner stability limits for e^ut = 0 (2 AU) and 0.25 (4 AU) 
are significantly different. This is despite the maximum inner bi¬ 
nary eccentricity being practically the same in the two cases. Qual¬ 
itatively, the difference in stability limits likely due to an increased 
eccentricity obtained by the planet, due to the eccentric third star, 
pushing the planet closer to the inner binary and causing instability. 

For fiout = 0.25 the planet remains stable for Op < 8 AU. For a 
circular tertiary the planet remains stable a little father out, before 
reaching instability at 9 AU. This is near where TKozai.p was seen to 
become shorter than Tprec.p in Fig. 5. These outer stability limits are 
on a significantly shorter-period than predicted by the circumpri- 
mary stability limit in Eq. 5 from Holman & Wiegert (1999): 35.54 
AU for gout = 0 and 24.54 AU for = 0.25. This is because Hol¬ 
man & Wiegert (1999) was calculated for a circular and coplanar 
planet, and not one that may potentially undergo Kozai cycles. 

In Fig. 7 we plot A/pp„ for the same simulations, and A/i„out 
in the bottom right figure. An ejection is often preceded by a large 
rise in the mutual inclination. For stable systems, there is a trend 
for the planets closer to the inner binary to remain close to copla¬ 
narity, whilst the planets farther away obtain at least a few degrees 
of mutual inclination. This means that even if Kozai cycles are sup¬ 
pressed, the tertiary star still has a perturbing effect. 



Figure 8 . Inner binary eccentricity over time for dp = 5 AU, eout =0 
and a varied Mp. Solid lines are from numerical integration and horizon¬ 
tal dashed lines are the eccentricity maxima calculated using Eq. 21. For 
nip = 0.0002Mo both the numerical and analytic results show ein.max = 0, 
i.e. a line across the bottom horizontal axis. 


5.2 What planets would suppress the KCTF process? 

The next question we have to consider is whether our binary orbit, 
with its primordial orbit of 100 d, will undergo KCTF without sup¬ 
pression by the planet. In Fig. 8 we plot over time for a planet 
at 5 AU with different planetary masses, in steps of 0.00005 Mq. In 
this simulation = 0. We compare these simulations with the an¬ 
alytic prediction from Eq. 21, showing that it is accurate to within 
roughly 5%. This small discrepancy may arise from Eq. 21 being 
derived under the assumption that = Oin Eq. 19. 

For a massless planet gin,max = 0.9345, and the expected final 
period of the inner binary is 5 d (Eq. 3). A planet of mass 0.00005 
Mq = 16.65 Mq causes a small reduction in gin.max = 0.9150. When 
nip is increased to 0.00015 Mq = 50 Mq the binary eccentricity is 
limited to 0.7875. In this scenario it is likely that the inner two stars 
never get close enough for tidal forces to have a noticeable effect 
over their lifetime, and hence there is likely to be no tidal shrink¬ 
age. A further increase in nip to 0.0002 Mq = 67 Mq completely 
suppresses the Kozai modulation. 

For planets at the inner (2 AU) and outer (8 AU) edges of 
the stability range, the critical planet masses for Kozai suppression 
from Eq. 22 are 5 Mq (0.016 Mj^p) and 300 Mq (0.94 Mj„p), respec¬ 
tively. 


5.3 A shrinking inner binary 

If the planet mass is sufficiently small such that the inner binary 
Kozai cycles are not suppressed, one would expect its orbit to 
shrink from ai„,i„it = 0.4827 AU (100 d) to ain,final = 0.0655 AU 
(5 d), according to the approximation in Eq. 3. As the inner bi¬ 
nary shrinks to 5 d the circumbinary precession timescale increases 
(Eq. 6). We demonstrate this in Fig. 5 by plotting Tpmc.p for multi¬ 
ple values of Pjn. The timescale turnover point shifts from Op = 7 
AU for a 100 d binary to rip = 3 AU for a 5 d binary. The Kozai 
timescale on the planet does not change as the binary shrinks, be¬ 
cause Pp and Pout are essentially static. 

We simulated all stable systems from Sect. 5.1 (flp between 2 
and 8 AU for Bout - 0 and flp between 4 and 7 AU for Bout = 0.25) 
with shorter binary periods at 25, 10 and 5 d, and looked for ejected 
planets. The simulations were over 25 Myr, which covers possible 
Kozai cycles for the planet but not for the shrunken inner binary, 
for which the Kozai timescale becomes very long. This integration 
time is considered reasonable because planets were already seen to 
survive high eccentricity excursions by the wider primordial binary. 
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Figure 9. Evolution of the planet’s eccentricity over 25 Myr around an inner 
binary of period 100, 25, 10 and 5 d, in the presence of an outer tertiary star 
at 100 AU. In the top figure Op = 4 AU and gout = 0. In the middle figure 
Op = 4 AU and gout = 0.25. In the bottom figure Up = 5 AU and gout = 0. 


Figure 10. Evolution of the mutual inclination between the planet and the 
inner binary, starting coplanar, over 25 Myr around an inner binary of period 
100, 25, 10 and 5 d, in the presence of an outer tertiary star at 100 AU, for 
the same simulations as in Eig. 9. In the top figure Up = 4 AU and gout = 0. 
In the middle figure Op = 4 AU and gout = 0.25. In the bottom figure Up = 5 
AU and gout = 0. 


which should have a more destabilising effect than shorter-period 
inner binaries. 

The three examples show: 

• For flp = 4 AU and gout = 0 the planet has a small eccentricity 
for Fiu = 100, 25 and 10 d. Around the primordial 100 d binary 
there is a bump in gp near 10 and 20 Myr, corresponding to Kozai 
cycles of the inner binary. This bump is not seen for shorter-period 
binaries because the Kozai cycles are longer than 25 Myr. When 
the binary has shrunk to 5 d the tertiary star’s strong influence on 


the planet, causes a significant raise in its eccentricity, yet it remains 
stable. There is an increase in the variation of A/p as decreases. 

• For flp = 4 AU and g^ut = 0.25 the results are qualitatively the 
same as in the case of a circular tertiary star. 

• For flp = 5 AU and gout = 0 the planet obtains a high eccentric¬ 
ity when the binary has reached a 10 d period, but remains stable. 
For a 5 d period binary the planet does not survive for more than 
a couple of million years, because the very short-period inner bi- 
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Figure 11. Numerically determined inner and outer stability limits for the 
planet as a function of Pin for eout = 0 (dark blue, solid) and 0.25 (light blue, 
dashed). The inner stability limit is defined from the test of the primordial 
binary in Sect. 5.1, and this does not change as the binary shrinks. 


nary is no-longer able to shield the planet from the perturbations 
from the tertiary star. This ejection is accompanied by a significant 
increase in A/p,in 

From the suite of simulations we derived rough stability lim¬ 
its as a function of Pin. In Fig. 11 we plot the inner and outer 
stability limits as functions of Pin, for eout = 0 (dark blue, solid) 
and 0.25 (light blue, dashed). Because the circumbinary precession 
timescale gets longer during the shrink, the tertiary star has an in¬ 
creasing influence over the planet. We found that this caused the 
outer stability limit to move inwards as Pin decreased. 

As Pin decreases the inner stability limit should also move in¬ 
wards (e.g. flcrit.cB = 0.14 AU for Pi,, = 5 d from Eq. 4). However, 
this is not a physically meaningful limit for any orbiting planets be¬ 
cause they had to form around the wider primordial binary where 
flcrit = 2 AU. Therefore, we define the inner stability limit as a con¬ 
stant, derived from the case of the primordial binary. 

Overall, the shrinking of the binary via KCTF can be seen as 
a destabilising process, as only planets formed relatively close to 
the inner binary have a change of surviving the shrinking process. 
There is only a narrow 2 AU region where a planet could possibly 
form and survive, assuming that the tertiary star has a circular orbit. 
For Sou, = 0.25 the inner and outer stability limits are the same 
(at 4 AU) when Pin reaches 5 d, and hence the likelihood of both 
forming a planet in these circumstances and having it survive the 
KCTF process down to a 5 d inner binary is small. 

6 ADDITIONAL CONSTRAINTS ON PLANET 
FORMATION AND EVOLUTION 

6.1 Protoplanetary disc environments 

Our analysis so far has been limited to n-body orbital dynamics. 
However planets are believed to form in discs and only under cer¬ 
tain favourable conditions. These necessary conditions further re¬ 
strict the possible range of disc radii that can allow planet formation 
in a stellar triple system. 

A circumbinary disc is expected to have a truncated inner edge 
near the n-body stability limit (Artymowicz & Lubow 1994). How¬ 
ever, it is considered theoretically challenging to form circumbi¬ 
nary planets close to this inner edge, due to secular forcing excreted 


by the binary that creates a hostile region for planetesimal accretion 
(Lines et al. 2014). The favoured theory is that planets are formed 
farther out in the disc in a more placid environment, before mi¬ 
grating inwards and halting near the inner truncation radius of the 
disc (Pierens & Nelson 2013; Kley & Haghighipour 2014). This 
theory naturally explains the pile-up of planets in Fig. lb. The fact 
that the abundance of circumbinary planets appears to be similar to 
that around single stars further supports this theory. A similar abun¬ 
dance implies a similar formation efficiency and environment, and 
a circumbinary disc only resembles a circumstellar disc far away 
from the inner binary. The implication is that a true inner limit for 
planet formation may be substantially farther out than the limit of 
Eq. 4. 

The formation of a circumprimary planet is also constrained 
by the interaction with the outer star, due to an outer truncation of 
the disc and perturbations on the planetesimals within. Like for the 
circumbinary case, the outer truncation radius in this case is also 
similar to the n-body stability limit (Artymowicz & Lubow 1994). 
However, here again the formation of planets in a circumprimary 
disc in the presence of a massive, inclined stellar companion can 
pose a significant theoretical challenge (Batygin et al. 2011). Theo¬ 
retical studies expect the particles in the protoplanetary disc to un¬ 
dergo significant modulations in eccentricity and inclination with 
different periods and amplitudes at different radial distances, mak¬ 
ing planetesimal collisions destructive rather than accretive (e.g. 
Marzari et al. 2009). This matches our results which show stable 
orbits, but with large fluctuations in eccentricity and inclination. In 
Fig. 12 we demonstrate this further with a zoomed version of the 
variation of A/p,in from Fig. 7 for Cp = 7 AU and Cout = 0.25. The 
fluctuations are of a couple of degrees in amplitude and on a short 
timescale with respect to planet formation. It is possible, however, 
that if these fluctuations are sufficiently mild they may be damped 
by the self gravity of the disc (Batygin et al. 201 1) or viscous dis¬ 
sipation (Martin et al. 2014). 

Like in the circumbinary case, the perturbations on the disc 
will be strongest near the outer truncation radius. This suggests 
a more restrictive outer formation boundary. These considerations 
are in accordance with observations showing that the properties of 
planets in circumprimary orbits are roughly indistinguishable to 
planets around single stars only if the binary companion is far¬ 
ther than ~ 100 AU (Duchene 2009), although the observational 
evidence for this is currently based on low number statistics. This 
finding is consistent with the scenario that giant planets are formed 
behind the snow line at at ~ 5 AU (e.g. Roberge & Kemp 2010), 
and that formation may be inhibited if the snow line is too close to 
or beyond the outer stability limit. 

6.2 Planet migration 

Our analysis so far has only considered planets with static 
semi-major axes. However, as previously discussed, the favoured 
paradigm is that the circumbinary planets generally migrate in¬ 
wards before being halted near the inner edge of the disc. The 
disc dispersal timescale (Alexander 2012) is expected to be much 
shorter than the KCTF timescale (Fabrycky & Tremaine 2007), and 
hence any migration will only occur around the primordial binary. 

This migration may have a positive or negative effect on planet 
survival, depending on the relative timescales of the inner binary 
Kozai cycles and disc migration. The disc truncation radius corre¬ 
sponds closely to the stability limit, which is in turn a function of 
fiin- If the planet migrates quickly and reaches the inner edge of the 
disc while the inner binary is still circular, then it will have migrated 
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Figure 12. Short-term variation in A/p in of a planet on an initially coplanar 
orbit at 7 AU around a primordial 100 d binary with a circular tertiary star 
at 100 AU, taken from the simulations in Fig. 7. 


in too far and will be ejected once ein is subsequently excited dur¬ 
ing Kozai cycles. Alternatively, if the planet migrates slowly then 
the inner binary will have already undergone a full Kozai cycle and 
the disc will be truncated farther out, meaning that the planet will 
not get too close. According to Sect. 5.3, having the planet near 
this edge will aid its later survival as the inner binary shrinks under 
KCTF. 


6.3 Multi-planetary systems 

Our analysis has assumed that there is only one circumbinary planet 
in the system, however with the discovery of Kepler-47 it is known 
that multi-planet circumbinary systems do exist in nature. Multi¬ 
planet stability has been analysed in circumbinary systems (Kratter 

6 Shannon 2013; Hinse et al. 2015) but never in the presence of a 
third star. 

When there are multiple planets in the system then there may 
be planet scattering events (e.g. Chatterjee et al. 2008 in the context 
of planets around single stars). This process may even be amplified 
in multi-stellar systems due to the eccentricity variations induced 
in the planets (like in Fig. 6 ), leading to more close encounters. 
Planet-planet scattering may lead to the orbits being pushed either 
inwards or outwards, potentially moving planets into unstable or¬ 
bits. A planet may also induce an apsidal precession on other plan¬ 
ets, which could possibly suppress Kozai perturbations from the 
tertiary star. 

If KCTF proceeds unhindered by the multi-planet system, the 
host binary will shrink and the secular timescale it imparts on the 
planets will sweep over a large range. This evolution would be sim¬ 
ilar to a sweeping secular resonance (Nagasawa et al. 2005): when 
the difference in binary-forced precession frequencies of the plan¬ 
ets matches the precession frequencies between the planets, it may 
lead to eccentricity excitation and perhaps destabilisation. 

7 SUMMARY AND DISCUSSION 
7.1 The general argument 

We here summarise our explanation for the lack of observed tran¬ 
siting planets around very short-period binaries, due to the KCTF 
process. Below is a list of all the key ingredients and constraints. 

(i) The planet forms around a long-period primordial binary un¬ 


dergoing high eccentricity Kozai cycles induced by the distant com¬ 
panion. 

(ii) To avoid ejection, the planet cannot be too close to the ec¬ 
centric inner binary, so it can only reside on a relatively wide orbit: 

^ ^^hin • 

(iii) However, the planet also cannot orbit too far from the inner 
binary, in order to have its Kozai cycles induced by the tertiary 
suppressed by the binary. 

(iv) If the tertiary star is eccentric, as is often the case, then the 
region of stability is even smaller. 

(v) If the planet is sufficiently massive and on a sufficiently close 
but stable orbit, it may inhibit the Kozai modulation of the binary, 
and therefore prevent KCTF shrinkage. 

(vi) Assuming that the binary orbit is able to shrink, some plan¬ 
ets become unstable, because the suppression of the Kozai cycle by 
the binary (constraint (iii)) becomes too weak. 

(vii) Even if Kozai cycles of the planet are suppressed by the 
inner binary, it still obtains some variation in A/p 

(viii) Considering the disc environment within which planets 
form leads to further restrictions being imposed, because the in¬ 
ner and outer edges of the truncated disc are probably too chaotic 
for planet formation. 

(ix) Finally, planet migration has the potential to either help or 
hinder the survival of planets. 

There is only a small parameter space where the planet can 
form, survive and not inhibit KCTF on the binary. In our example 
in this section only roughly Neptune-mass planets between 2 and 
3 AU fulfilled all of the above criteria. Furthermore, the surviving 
planets were misaligned by over 10°. This triple star configuration 
would be even more restrictive if the perturbing strength of the ter¬ 
tiary were increased by (i) an increased mass ratio M 3 /(Mi Mp), 
(ii) a decreased semi-major axis ratio flout/oin: (fit) an increased 
or (iv) an increased mutual inclination A/i„ out. leading to higher ec¬ 
centricity Kozai cycles. 

We conclude that most triple star systems evolving under 
KCTF are not conducive to hosting planets. Alternatively they host 
planets biased towards small masses, long periods and misaligned 
orbits, which are difficult to detect via transits. 

The known circumbinary planets around wider binaries likely 
formed in a more placid environment, suggesting that companion 
stars either do not exist or are too far away to have an effect, like in 
PH-l/Kepler-64. We await future observations to confirm this. The 
fact that these planets have been found near the inner stability limit 
suggests that the binary orbit has not shrunk over a long timescale 
(e.g. via KCTF) after planet migration has finished. 

Our analysis could be further extended by testing a wider 
range of orbital parameters to quantify different regimes of secular 
evolution of planets in triple star systems, similar to the approach 
taken in Takeda et al. (2008), but this is beyond the scope of the 
current more qualitative investigation. 

Whilst this general argument developed may account for the 
majority of the non-detections of transiting planets around short- 
period binaries, it may not be the only effect present. Very close bi¬ 
naries are tidally locked, which increases the rotation speed and can 
lead to increased stellar activity. The standard Kepler 30-minute ca¬ 
dence may lead to insufficient sampling at the shortest periods. And 
finally, Martin & Triaud (2015) calculated that some of the closest 
eclipsing binaries may be sufficiently inclined with respect to our 
line of sight that transits by coplanar planets are not geometrically 
possible. 
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7.2 Looking ahead 

Some level of planet formation and survival might be possible in 
this shrinking binary scenario, however the bias towards planets 
with small masses and wide, misaligned orbits creates significant 
detection limitations. To find such planets it is likely that new tech¬ 
niques will need to be developed. The four years of Kepler pho¬ 
tometric data are yet to be exhaustively searched, and may still 
yield discoveries of misaligned planets via eclipse timing variations 
(Borkovits et al. 2011), or transits on non-eclipsing binaries (Mar¬ 
tin & Triaud 2014), particularly on binaries found with the BEER 
technique (Faigler et al. 2012). A re-observation of the Kepler field 
by the future PLATO satellite may lead to transits on the existing 
sample of eclipsing binaries by new planets that have precessed into 
view. Radial velocity surveys can detect misaligned circumbinary 
planets, however short-period binaries are tidally locked, leading to 
an increased rotation velocity that decreases spectral precision. 

There is also potential for gravitational lensing and direct 
imaging (Delorme et al. 2013; Thalmann et al. 2014) discoveries 
of long-period circumbinary planets. Finally, GAIA astrometry will 
be sensitive to potentially hundreds of giant circumbinary planets 
on periods of a few years, and can even provide a direct measure¬ 
ment of the mutual inclination (Sahlmann et al. 2014). 

Continued observations, whether they reinforce this dearth or 
lead to surprising new discoveries, will allow this relatively new 
problem in exoplanetary astrophysics to shed new light on a funda¬ 
mental field of stellar physics. 
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